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Abstract 

We study CP violation in the lepton sector of the supersymmetric extension of the Stan- 
dard Model with three generations of massive singlet neutrinos with Yukawa couplings Y v 
to lepton doublets, in a minimal seesaw model for light neutrino masses and mixing. This 
model contains six physical CP- violating parameters, namely the phase S observable in oscil- 
lations between light neutrino species, two Majorana phases 0i j2 that affect (3(3q v decays, and 
three independent phases appearing in Y„Yj, that control the rate of leptogenesis. Renor- 
malization of the soft supersymmetry-breaking parameters induces observable CP violation 
at low energies, including T-odd asymmetries in polarized /i — > eee and r — > III decays, as 
well as lepton electric dipole moments. In the leading-logarithmic approximation in which 
the massive singlet neutrinos are treated as degenerate, these low-energy observables are 
sensitive via YjY v to just one combination of the leptogenesis and light-neutrino phases. 
We present numerical results for the T-odd asymmetry in polarized fi — > eee decay, which 
may be accessible to experiment, but the lepton electric dipole moments are very small in 
this approximation. To the extent that the massive singlet neutrinos are not degenerate, 
low-energy observables become sensitive also to two other combinations of leptogenesis and 
light-neutrino phases, in this minimal supersymmetric seesaw model. 
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1 Introduction 



The solar |T], and atmospheric || neutrino anomalies, which imply the existence of non- 
zero masses for the light neutrinos, provide the first experimental evidence for the existence 
of physics beyond the Standard Model (SM). A minimal extension of the SM includes three 
very heavy singlet neutrinos iV?, whose Yukawa couplings Y v to the light neutrinos explain 
naturally the smallness of their masses, via the seesaw mechanism ||. At the same time, 
the electroweak scale must be stabilized against large radiative corrections. In particular, 
after introducing right-handed neutrinos, a quadratically-divergent contribution to the Higgs 
boson mass proportional to has to be cancelled. This is most commonly achieved 
by supersymmetrizing the theory, leading to the minimal supersymmetric extension of the 
Standard Model (MSSM) with singlet neutrinos. 

Neutrino-flavour mixing originates from off-diagonal components in the Yukawa interac- 
tion N C Y V LH2, in a basis where the charged-lepton and singlet-neutrino mass matrices are 
real and diagonal. Renormalization effects due to this interaction also induce flavour mix- 
ings in the soft supersymmetry-breaking slepton mass terms |J. This may lead to observable 
rates for charged-lepton flavour- violating (LFV) processes such as /i — > cy, fi-e conversion 
in nuclei, /i — > eee and r — > 3£ [7L |SL Bj, where £ — e, fM denotes a generic light charged 
lepton. LFV is also observable in principle in rare kaon decays, but at rates that are likely 
to be far below the current bounds jlQ| . 

In general, Y v is complex, leading to CP violation in neutrino oscillations and in the 
induced rare LFV processes, as well as in Majorana phases for the light neutrinos and in 
electric dipole moments for the charged leptons. The existence of CP violation in Y v is also 



required if the observed baryon asymmetry in the Universe originated in leptogenesis [11 



The purpose of this paper is to clarify the relations between these different manifestations 
of CP violation in the lepton sector, and to present numerical estimates of the T-odd CP- 
violating asymmetry At in /i — > eee decay, the electric dipole moments of the electron and 
muon. We argue that measurements of CP violation using charged leptons, combined with 
CP violation in the light-neutrino sector, in principle enable the leptogenesis phases to be 
extracted - within the framework of the minimal supersymmetric seesaw model. 

If the solar-neutrino mass-squared difference Am^ and the element U e z of the Maki- 
Nakagawa-Sakita (MNS) neutrino-mixing matrix U are not too small, the CP- violating phase 
S in U, which is analogous to the Cabbibo-Kobayashi-Maskawa (CKM) phase in the quark 
sector, can be measured via CP- and T-violating [12j observables in neutrino oscillations 
using a neutrino factory or possibly a low-energy neutrino superbeam. The recent SNO 
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result |2| encourages this possibility, since it further favours the large-mixing-angle (LMA) 
solution to the solar- neutrino deficit flIB| . 

As mentioned above, processes that violate charged-lepton flavour can provide important 
complementary information on the leptonic CP-violating phases. These may be measured 
using intense sources of stopped muons. The SINDRUM II experiment is designed to be 



sensitive to B([iTi —* eTi) ~ 10 14 [14], and the MECO project would be sensitive to 



-14 



is 



B(fiAl — > eAl) ~ 10~ 16 fig] . The experiment with the sesitivity Br(fi — > e^y) ~ 10 
proposed at PSI The PRISM project || [17| and the front ends of neutrino factories now 
under consideration at CERN [18] and elsewhere will provide beams of low-energy muons 
that are more intense by several orders of magnitude than any of the present facilities. This 
will enable the construction of stopped-muon experiments able to probe LFV processes with 
sensitivities Br(fi — > cy) ~ 10 -15 , Br{^i — > eee) ~ 10~ 16 . The latter sensitivity opens the 
way to measuring the T-odd, CP- violating asymmetry At{h —>■ eee). 

A measurement of the CP-violating electric dipole moment (EDM) of the muon with a 
sensitivity <i M ~ 5 x 10~ 26 e cm would also be possible [18|. However, because the Yukawa 
coupling constants Y v appear in the renormalization-group equations (RGEs) only in the 
Hermitian combination YjY u , CP-violating phases are induced only in the off-diagonal terms 
of the slepton masses. This implies suppression of the EDMs of the electron and muon, 
whereas CP violation may occur in full strength in charged LFV processes, such as fi — > eee [J 

Another arena to probe LFV is provided by rare r decays. There has been some discussion 
in the literature of r — > £7 decays [P7[ , and one could in principle hope to measure CP- 
violating asymmetries in the various r — ► 3£ decays. Another possibility is to search for LFV 



in sparticle decays p8[ |2~9fl , e.g., Xi ~ > e f i Xp where CP-violating asymmetries analogous 
to At(h —>■ eee) can also be defined in principle. However, we do not investigate these 
possibilities further in this paper. 

We concentrate here on CP-violating observables in the /i sector, assuming that the 
only sources of LFV and CP violation are the interactions with heavy singlet iV c neutrinos. 
We start by discussing general parametrisations of the Yukawa matrix Y v in terms of the 
high- as well as the low-energy observables, paying particular attention to the counting of 
physical degrees of freedom and their relations to CP-violating observables. Subsequently, we 
analyse the renormalization-group running of soft supersymmetry-breaking terms, assuming 
universal boundary conditions at the GUT scale. Our first objective in this analysis is to 



1 In the light of very stringent constraints from electron, neutron and mercury EDMs |19 



2§ we 

neglect the possible phases in diagonal soft supersymmetry-breaking terms throughout this paper. This is 
natural in mechanisms which generate only real soft terms, such as gravity- p2|, gauge- [^3|, anomaly- [ p4[ . 
gaugino- |2q| and radion-mediation pq| mechanisms. 
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demonstrate in principle the complementarity of the different observables, to see how all 
the CP-violating phases of the minimal seesaw model come into play, and to clarify the 
relationship of the observable phases to the phases appearing in leptogenesis. In the leading- 
logarithmic approximation, in which the heavy singlet neutrinos are treated as degenerate, 
this renormalization is sensitive to just one combination of the leptogenesis and light-neutrino 
phases, but two other combinations contribute beyond this approximation. We illustrate our 
results in a simple two-generation model. We then present numerical estimates of Br{ji — > 
ej), Br(/i — > eee) and A T (/i — > eee), taking into account the present knowledge of neutrino 
mixings and masses as well as bounds on sparticle masses. We find that the magnitude 
of At is in general anti-correlated with the rate of /i — > ej, and may be large in some 
models compatible with the experimental upper limit on \i — > e^ decay. If a cancellation 
occurs between different contributions to the ji — e — 7 vertex, so that the box and penguin 
diagrams contributing to /i — > eee become comparable in magnitude, the T-odd asymmetry 
At in fj, — > eee may be as large as ~ 10%, while Br{ji — > eee) remains appreciable. However, 
the EDMs of the \i and e are rather small in the minimal seesaw model. 

It is important to note that the neutrino-oscillation phase 5 and the Majorana phases 
(pi t 2 are completely independent of the three physical phases in the quantity Y u Yj that enters 
in leptogenesis calculations. On the other hand, A T and the other renormalization-induced 
observables depend on mixtures of the light-neutrino and leptogenesis phases. Thus, neutrino 
factories and LFV measurements provide complementary information on the leptonic CP- 
violating phases. In particular, observation of At is possible even if CP violation in neutrino 
oscillations is unobservable, i.e., if either 5 = 0, U e3 = or the solar-neutrino deficit is not 
explained by the LMA solution. However, in the minimal supersymmetric seesaw model, it is 
possible that a combination of CP-violating observables in the neutrino and charged- lept on 
sectors may provide constraints on the angles and phases responsible for leptogenesis. 

Our work is organized as follows. In Section 2, we consider general parameterisations 
of the neutrino Yukawa couplings Y u and discuss CP violation in the minimal supersym- 
metric seesaw model. In Section 3, we give general formulae for the EDMs of the charged 
leptons, ji — > e7, and ji — > eee, including the latter's T-odd asymmetry A T . We present the 
results of the numerical analysis in Section 4. Finally, Section 5 is devoted to a discussion 
and our conclusions concerning the observability of the CP-violating phases in the minimal 
supersymmetric seesaw model. 
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2 CP Violation in the Lepton Sector of the Minimal 
Supersymmetric Seesaw Model 



We consider the MSSM with three additional heavy singlet-neutrino superfields N c i, consti- 
tuting the minimal supersymmetric seesaw model. The relevant leptonic part of its super- 
potential is 

W = N^(Y l/ ) ij LjH 2 - E?(Y e )ijLjHi + ^N^M^N* + ^H 2 H X , (1) 

where the indices i, j run over three generations and A4ij is the heavy singlet-neutrino mass 
matrix. Taking account of the possible field redefinitions, this minimal supersymmetric 
seesaw model contains 21 parameters: 3 charged-lepton masses rri£, 3 light-neutrino masses 
M.® , 3 heavy Majorana neutrino masses M. D , 3 light-neutrino mixing angles 0y : 1 < 
i 7^ j < 3, 3 CP-violating light-neutrino mixing phases 5, 0x,2 (the MNS phase and two 
Majorana phases), and 3 additional mixing angles and 3 more phases associated with the 
heavy-neutrino sector. 

2.1 High-Energy Parametrization 

In order to clarify the appearance and roles of these parameters, we first analyze ([!]) in a 
basis where the charged leptons and the heavy neutrinos both have real and diagonal mass 
matrices: 



{Y e ) l3 = Y°8„ 



My = MfSn, (2) 

where Ai D = diag(Mjv 1 , M^ 2 , Mat 3 )- A "priori, the neutrino Yukawa- coupling matrix Y v has 
nine phases, which can be exposed by writing it in the form: Y v = Z*Y^X^ P* , where Y® is 
diagonal and P\ = diag(e* CTl , e" 72 , e" 73 ). However, in the basis (0) one may redefine the left- 
handed lepton fields L i; and thus rotate away the three phases in Pi, which are unphysical. 
Thus the Yukawa-coupling matrix may be written in the form: 

(Xu)ij = Z* k Yj^X kj . (3) 

The matrix X is the analogue in the lepton sector of the quark CKM matrix, and thus it 
has only one physical phase. On the other hand, we can always write Z in the form 

Z = P{ZP 2 , (4) 
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where Z is a CKM-type matrix with three real mixing angles and one physical phase, and 
Pi ; 2 = diag(e jei ' 3 , e^ 2 ' 4 , 1) are diagonal matrices containing two phases each. Thus Z has 5 
physical phases to add to that in X, and all six real mixing angles and six phase parameters 
in this basis are physical observables. 

We now study the combination Y V YJ of the Yukawa couplings, which governs leptogenesis 
in this minimal seesaw model. It is straightforward to see from (|3]) that 

Y V Y} = PfT(Y v D ) 2 Z T P 1 , (5) 

which depends on just three of the CP-violating phases, namely the two phases B 12 in P\ 
and the single residual phase in Z, as well as the three real mixing angles in Z . This is 
consistent with the observation that, since the overall lepton number involves a sum over 
the light-lepton species (both charged leptons and light neutrinos), one would not expect 
leptogenesis to depend on the 6 MNS angles and phases. 

On the other hand, as we discuss in more detail below, mixing and CP violation in the 
slepton sector of this minimal supersymmetric seesaw model is controlled by the combination 
YjY u of the neutrino Yukawa couplings, in the leading-logarithmic approximation where 
Mqut 3> Mat 12 3 . It is again straightforward to see from (|3|) that 

YjY u = X{Y u D fxK (6) 

Therefore, in this approximation, CP violation in charged LFV processes arises only from 
the one physical phase in the diagonalizing matrix X. 

2.2 Low-Energy Parametrization 

We now reconsider leptonic CP violation from a more familiar point of view fl30| , namely 
that of the effective low-energy theory obtained after the heavy neutrinos are decoupled. In 
this energy range, physics is described by the following effective superpotential: 

W eff = UH 2 (yJ (M D Y l Y v ) L ] H 2 -E<r(Y e ) i] L J H 1 , (7) 
V / ij 

where the effective light-neutrino masses are given in the basis (0) by 

M u = Yj (M D ) ^ Y u v 2 sin 2 p, (8) 

where v = 174 GeV and as usual tan/3 = v-ijv\. The mass matrix M. v can be diagonalized 
by a unitary matrix U : 

U T M U U = M° , (9) 
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where M.® = diag(m^, m^, m„ 3 ). Since M. v is a symmetric matrix and contains in general 
six phases, U must also have 6 phases. It can be expressed in the form 

U = P 2 VP , (10) 

where P$ = diag(e~ 1 ^ 1 , e~ 1 ^ 2 , 1), P 2 = diag(e lcn , e za2 , e"* 3 ) and V is the MNS matrix written 
in the CKM form: 



/ C13C12 C13S12 s 13 e lS 



V 



-C 23 S 12 - S 23 S 13 C 12 e t5 C23C12 - S 2 3Sl3Si 2 e i<5 S23C13 J- (11) 
S 23 S\ 2 — C 23 Si 3 Ci 2 e %S — S23 c 12 — C23Si3Si2e i<5 C23C13 / 

The phases in P 2 (JD) can be removed by redefinition of the Li fields, leading to a new basis 
in which 

U = VP . (12) 

This differs from the basis (Q) by the phase rotation P 2 . The new basis is appropriate if 
one works with the effective low-energy observables in the effective superpotential (|7|), e.g., 
for studying neutrino oscillations. Indeed, the mixing angles : 1 < % j < 3, whose 
sin, cos we denote by Sy, c^-, are measurable in neutrino-oscillation experiments, as is the 
CP-violating MNS phase 6. One combination of two CP-violating Majorana phases 0i 5 2 is 
in principle measurable in /3/3o u experiments. 

The physical interpretation of the Yukawa couplings in (|l|) is made more transparent in 
the basis @, which does not contain the unphysical low-energy phases in P 2 that we rotated 
away in the previous paragraph. Note that one must change X — > P 2 X if one works in the 
basis (HI). 

Our objective in this paper is to study CP-violating observables which are sensitive to 
different physical phases. For this purpose, we need a proper parametrization of the input 
parameters of the model. The most straightforward choice is to work in the basis (0) and 
to choose the physical observables in (0) and (0) as the input parameters. In this case, the 
physics is entirely transparent. However, the present experiments do not measure heavy- 
neutrino masses, their Yukawa- coupling and mixings directly. All the information we have 
on neutrinos comes from the low-energy neutrino-oscillation and /3/3q u experiments. If we 
choose the input parameters from (^|), (||]), we have to check every time that the induced 
M. v in (H) agrees with the experimental data. Instead, one can attempt to use the effective 
low-energy observables as an input. 

To this end, we first rewrite the seesaw mechanism in the different form: 

R = VM^^Y v uJm^' 1 : R T R=1, (13) 
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which is equivalent to (|]). Starting with any given Y v and Ai D as input parameters, we 
obtain as outputs the seesaw- induced low-energy parameters Ai® and U, and an auxiliary 
complex orthogonal matrix R. It is possible to choose different parameter sets for Y v and 
M D that give the same low energy effective M.® and U, but lead to different values for R. 

One can turn the argument around |3J, and parameterize the neutrino Yukawa-coupling 
matrix in terms of an arbitrary complex orthogonal matrix R' as follows: 



Yl = (14) 

vsm p 



We emphasize that the output Y' v in this parametrization is in the low-energy basis (12|), and 
therefore contains unphysical phases. If one wants to use the induced Yl to parametrize the 
superpotential (j^), one should be careful to count correctly the physical degrees of freedom. 
We now form the combinations Y'JY'J and Y'^Yl out of (|Hp. In the first case, we obtain 



which contains three independent physical phases that are given entirely in terms of the pa- 
rameters in the orthogonal matrix R' . This is consistent with (|5|), and the new parametriza- 
tion therefore has not changed the counting of phases in Y u Yj. On the other hand, we also 
obtain from ([[¥) 



Y'JYI = 2 . V " yt. (16) 

v z sm p 

This expression also appears to contain three phases, which are combinations of all the 
parameters in U and R' . 

However, according to (|J), Y u Yj is supposed to contain only one physical phase. What 
has happened? The answer is that physics has not changed, and thus two out of three phases 
in YjYl are unphysical. This is the case because we are working in the low-energy basis 
(|12|), and not in the basis (Q). The three phases in P2, which were rotated away in defining 
U, appear now in Y' v . Instead of (||), we now have Y'jYl = P 2 X(Y V D ) 2 X^P*. One overall 
phase is irrelevant, and the two unphysical relative phases in P 2 explain the faulty phase 
counting in fllSp. 

In the following, we show explicitly that the unphysical phases in P 2 cancel out in the Jarl- 
skog invariants which can be constructed using YjYl. Therefore, in the leading-logarithmic 
approximation, all the CP-violating LFV observables depend only on the one physical phase 
in (|16|) , which is a combination of the phases in U and R'. Henceforward, we omit the 
superscript ' , but one must still be careful to distinguish between the different bases. 
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2.3 Relations to CP-Violating Observables 

So far we have only considered the parametrization of the input neutrino parameters, which in 
general are complex, and the 6 resulting independent CP- violating phases. We now consider 
how physical observables depend on these various phases. 

2.3.1 Leptogenesis 

At present, our only experimental knowledge on CP violation in the lepton sector may be 
obtained from the baryon asymmetry of the Universe, assuming that this originated from 
leptogenesis. In leptogenesis scenarios, initial B — L asymmetries e l appeared in decays of the 
heavy neutrinos N? in the early Universe, as results of interferences between the tree-level 
and one-loop amplitudes for N% decays. The L asymmetry in the decay of an individual 



species Nf is given in the supersymmetric case [31] by 



-E 

8«r 



Im 






li - 




1 V 


2 



V£T[Log(l + l/xi) + 1——^} , (17) 



where x\ = (MnJM^) 2 and both triangular and self-energy type loop diagrams are taken 
into account. This L asymmetry is converted into the observed baryon asymmetry by 
sphalerons acting before the electroweak phase transition. It is clear from (|T7| ) that the gener- 
ated asymmetry depends only on the phases in Y^Yj. Hence, according to the parametriza- 
tion (0), the only phases entering in the calculation of the baryon asymmetry of the Universe 
are those in R. In order to demonstrate the feasibility of leptogenesis, it would be necessary 
to prove that at least one of the phases in R is non-zero. Moreover, as we shall see, at least 
one of the real part of the mixing angles in R must also be non-zero, and one would need to 
control other parameters, such as the heavy-neutrino mass spectrum, before being able to 
calculate the baryon asymmetry in terms of R, or vice versa. 

2.3.2 CP Violation in Neutrino Oscillations 

Measuring this is one of the main motivations for building neutrino factories. We assume 
that the real MNS mixing angles 0^ and the mass-squared differences 5m?- = m^. — m 2 . are 
all non-vanishing, in which case the the MNS phase S in (|TT[) is in principle observable. It 
is, realistically, observable in long-baseline neutrino factory experiments if the LMA solution 
of the solar neutrino problem is correct. The Majorana phases <f>\p do not affect neutrino 
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oscillations at observable energies, but do affect (3(3q v decay. The conventional nuclear (3(3q v 
experiments measure one combination of the light-neutrino masses m Vi and the Majorana 
phases 01^- As in the CKM case, one can introduce a Jarlskog invariant that characterizes 
the strength of CP violation in neutrino oscillations: 

= 5m\ 2 Sm 2 23 6m 2 31 Im [V X1 V^V 22 V 2 \] . (18) 

One sees explicitly that the Majorana phases (f>i^ cancel out in J u . 

It is clear from (|8|) that, from the high-energy point of view, 5 depends on all the six 
independent phases in Y u , including those in the combinations Y u Yj and YJY U . On the other 



hand, in the low-energy parametrization of (fH]), the phase 5 is taken as an input parameter 



2.3.3 Renormalization of Soft Supersymmetry-Breaking Terms: 
Flavor-Changing Processes 

In the minimal supersymmetric seesaw model, renormalization induces sensitivity to the 
neutrino Yukawa couplings Y v in the soft supersymmetry-breaking parameters in the slepton 
sector, in particular to the CP- violating phases in Y v . These may have measurable effects 
on several CP-violating lepton observables, including asymmetries in LFV decays, which 
are observable in rare fi and/or r decays, and electric dipole moments. The fi electric 
dipole moment as well as rare fi decays may be measurable using slow or stopped muons 
produced at the front end of a neutrino factory. In this subsection, we concentrate on the 
flavor- changing processes, such as asymmetries in LFV decays, and we will discuss flavor- 
conserving processes, such as the electric dipole moment. 

The soft supersymmetry-breaking terms in the leptonic sector of the minimal supersym- 
metric seesaw model are 

+ (N^A^L^ - Et{A e ) l3 l 3 E x + l -Nr{B N ) lJ N'r 

+-M X BB + l -M 2 W a W a + ~ M 39 a 9 a + h.c^j . (19) 

We assume that the soft supersymmetry-breaking terms have universal boundary conditions 
at the GUT scale M GUT ~2x 10 16 GeV: 

( m l)ij = ( m l)ij = ( m %)ij = mil, 
m \ x = m H 2 = m o , 
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(A e )ij — Ao(Y e )ij , [A. v )ij — Ao(Y u )ij , 

M 1 = M 2 = M 3 = mi/a . (20) 

At lower energies below Mqut and above the heavy-neutrino mass scale Mat, which we 
assume to be <C Mqut, off-diagonal entries in Y v generate via the renormalization-group 
running off-diagonal entries in the effective soft supersymmetry-breaking terms. In the 
leading-logarithmic approximation the flavor- dependent parts of the soft supersymmetry- 
breaking terms are given by 



GUT 



(Sm%) ~ -^ 2 (3rnl + Al)(Y e Y e %log MGUT 



» 4vr 2V u u/v c ejJ0 ° M N ' 

(M e k « -^^(SF^ + Fj^^log^f^. (21) 



Here, the Yukawa coupling constants are given at Mjy, and then Y" e is diagonal. This means 
that ml remains diagonal in this approximation. Below Mat, the heavy neutrinos decouple, 
and the renormalization-group running is given entirely in terms of the MSSM particles 
and couplings, and is independent of Y v . We use in our numerical examples full numerical 
solutions to the one-loop renormalization-group equations, but the approximate analytical 
solutions ([H]) are useful for a qualitative analysis. 

It is important to notice that, in the leading-logarithmic approximation (pi]), the only 
combination of neutrino Yukawa couplings entering the renormalization-group equations is 
YJY U . This implies that CP-violating phases are induced only in the off-diagonal elements of 
(m|)jj and (A e )ij, and further indicates that the lepton-flavour conserving but CP- violating 
observables like the electric dipole moments of charged leptons are naturally suppressed [p0 |, 



while CP violation in the charged LFV processes should occur in full strength. This is 
analogous to CP violation in the quark sector of the Standard Model, which is also directly 
related to flavour-changing processes. As we saw earlier (|GJ), the combination YJYv depends 
on just one CP-violating phase, namely that in the matrix X. Therefore, in the leading- 
logarithmic approximation, all slepton-induced observables are independent of the phases 
associated with leptogenesis, which are combinations of those in the matrices Z and P\ (|B|), 
in the high energy parametrization. On the other hand, in the low-energy parametrization, 
YJYv depends on one combination of the phases in U and R, as explained in subsection 2.2. 

Since YJYv depends on only one physical phase, there is only one invariant for m| de- 
scribing the strength of CP violation in any process induced by sleptons. By analogy with 
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the Standard Model quark sector, this can be taken to be 



Jf 



Im 



H) 12 H) 23 K) 



31 



Additional invariants including the A terms can be constructed: 



Im 



W 12 H) 23 (m|) 3i 



(22) 



(23) 



and cyclic permutations, and similar invariants with two or three (^4 e )y factors. However, 
in this model they are all related to the basic invariant (p2|) , and proportional to 



Im 



31 



6Y U \ 2 6Y U \ 3 5Y U 2 31 lm[X u X* 2 X 22 X^] 



(24) 



in the leading- logarithmic approximation ( pi]) . Here, 5Y^- = (}^P) 2 — (V^. 

The above analysis is modified when one includes in the renormalizat ion-group running 
effects associated with the non- degeneracy of the heavy neutrinos: ^ ■ In this case, 
{8rrv~)ij in (|21| ) is replaced as follows: (<5m|)jj — > (5m|)y + (<5m~J .., and u5>m|J . is given 
by 



(5ml) 



M, 



:(3m 2 + A 2 )(Y}LY„) tJ . : L = log -^A 



J7r2 v--u ■ "o/v-v—^v °M^. iJ ' 



(25) 



where is now interpreted as the geometric mean of the heavy singlet-neutrino mass 
eigenvalues M^ t . The first term in fl25|) contains the matrix factor 



y ] ly = xy d p 2 z t lz*p;y d x\ 



(26) 



which induces some dependence on phases in ZP 2 . In the three-generation case, there are 
two independent entries in the traceless diagonal matrix L, so the renormalization induces 
in principle dependences on two new combinations of these phases, as well as the single 
phase in YjY„. Thus low-energy observables become sensitive to all three leptogenesis 
phases. However, the dependences on the two extra phases are suppressed to the extent 



that log ^« log 



M, 



M N . 



OUT 



M N 



2.3.4 Renormalization of Soft Supersymmetry-Breaking Terms: 
Flavor-Conserving Processes 

As mentioned above, since the CP-violating phases are in the off-diagonal components of 
the soft supersymmetry-breaking terms, the electric dipole moment of lepton is naturally 
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suppressed. The following is the lowest-order combination of the Yukawa couplings Y v and 
Y e whose diagonal components have imaginary parts: 



J, 







edm 



Im 



Y P Y},Y„ [yjv , vtvl vtv 



1/ 1 V 



■ 1/ - 1 f 



2Y ei 6Y e % 6Y U \ 2 5Y? 23 6Y? 31 Im^^X*]^, 



■ijk 



(27) 



and the dominant contributions to the electric dipole moment are proportional to it0. Since 
YjY u in J c ^ m comes from the radiative correction to the soft supersymmetry-breaking terms, 
the leading contribution to the electric dipole moment is proportional to log 3 M ^ I v N r when 
log < 4tt. The dependence on Y}Y e in J e ( J m comes from the radiative correction in the 
soft supersymmetry-breaking terms or the tree-level mass matrix of the charged sleptons. In 
this subsection, we present the Jarlskog invariant for the soft supersymmetry-breaking terms 
contributing to the electric dipole moment in (9 (log 3 M ^ [ VT )- Also, we discuss cases where 



this approximation is invalid, namely when i) tan/? 3> 1, or ii) non-degeneracy between the 
heavy singlet neutrinos induces dependences of m| and A e on phases in the product ZP 2 . 

In order to evaluate the contribution to the electric dipole moment in 0(\og 
need the corrections to the soft supersymmetry-breaking terms at O (log 2 M ^ T , 



3 M. 



GUT 



we 



which are 



(4tt) 



M, 



GUT 



M, 



N 



2) 2 
>m E 



{6® A 



ejij 



(47T) 

0. 



; A 2 (3Y e YXY^ + Y e YX„Y}) 3l log 1 



M, 



GUT 



M, 



N 



(28) 



Here, we neglect irrelevant terms with a trace over flavor indices, or which are flavor- 
independent. The Yukawa couplings are evaluated at M N . From these equations and (|2~ID , 
non- vanishing contributions to J® m arise from the following combinations of O (log 3 M F f UT ) : 



Im 



8A P 5A\5A f 



rA 3 jillog 3 



Mi 



GUT 



Im 



5A P YlYJ {2) rn 



(47r) 6 dm 

12 4 3 7 (i) W 

1 J odm 10 S 



(4vr; 



M N 
Mqut 
M N ' 



(29) 
(30) 



In ([II]), the combination YJYe arises from the tree-level mass matrix of the charged sleptons. 
It is found from (^,^) that the electric dipole moments depend strongly on A Q and less on 
m . 



2 Similar studies for the electric dipole moment of neutron in the MSSM are done in [B2|, assuming that 
all CP violating phases come from the CKM matrix. 
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When tan/5 ^> 1, ( |30]) is proportional to tan 2 /3 and ( p9|) is not enhanced. On the other 
hand, terms such as 



Im 



(31) 



are proportional to tan f3 and log £j? T . Thus, they may make sizeable contributions to 
the electric dipole moments for tan (3 1, even if log M M v N r ^ 47r. 

If the heavy neutrinos are not degenerate in mass, they induce dependences of the soft 
supersymmetry-breaking terms on phases in Z/P2, as mentioned in the previous Section, 
which then contribute to the electric dipole moments. The Jarlskog invariant J^ m depends 
on Y^Y e , and this factor suppresses the electric dipole moment when tan/3 is small. In this 
case, the non-degeneracy of the heavy neutrinos may have a more important effect on the 
electric dipole moment. The corrections of 0(\og M M ^ T log j^r) are 

GUT 



~ J^Ai A ° Y ^ n W LY "> Y ^ + 7 ^ LY -> Y ' Y "\k log ^Jjf ■ (32) 

Here, we neglect terms with Yj>Y e factors. The interesting point is that the second term 
111 can have imaginary parts in the diagonal components, and thus can contribute to 

the electric dipole moment Since phases in S^A e arise from ZP 2 , we do not need three 
generations of leptons in order for 8^A e to have imaginary parts in the diagonal terms. The 
behaviour of this contribution will be discussed in the next subsection. 

2.4 Two- Generation Model 

We now demonstrate the interdependences of the above physical observables in a toy two- 
generation model. In this model, X has no physical phase while there may be one phase in 
Z. We parametrize the light- and heavy-neutrino masses and R as follows: 




(33) 



3 Whilst the combination (6A e 5m~)ii has an imaginary part, it does not contribute to the electric dipole 



L 

moments, since Im[(<L4 e + 5A e )(5m 2 ^ + 8rrv~J\a = 0. 
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In this model, the leptogenesis invariant is Im {/Y v Yi/^ (Y v Yj^ 

r / i\ 12 n 



2 xlm 



12 



X 



Re 



(Y„Yj) 



where 
Im 
Re 



Y Y 1 

1 V 1 V 



(Y„Yj) 



12 



12 



2v 2 sin 2 /3 
{m Ul -m V2 )yfMjA 2 



so that 



Im 



(Y y Yj) 21 (Y u Yj 



21 



2v 2 sin 2 (3 
{ml-ml 2 )M 1 M 2 



sin26L 



sinh26',sin2^ r 



(34) 



(35) 



2w 4 sin 4 /3 

As explained above, the phase 0, in R controls leptogenesis, and the mixing angle 9 r must 
also be non-vanishing. 

As concerns neutrino observables, we recall that there is no analogue of the MNS phase 
5 in this two-generation model. There is one CP-violating Majorana phase <fi for the light 
neutrinos, but this does not contribute to leptogenesis, as we argued previously on general 
grounds and now see explicitly in (|35|) . 

We now consider the quantity YjY„ which controls the renormalization of the soft super- 
symmetry-breaking terms in the leading-logarithmic approximation, in particular (YjY u ) 12 
which has a non-zero imaginary part. For illustrational purposes, we assume that the light- 
neutrino mass matrix has maximal mixing: 



U 



l/y/2 1/V2 \ 
-1/V2 l/y/2 J 



1 



(36) 



where is a light-neutrino Majorana phase. In this case, 



Re 



Im 



(YjY u ) 



Y T F 

1 V 1 V 



12 



m V2 + m 



12 



4v 2 sin P 2 

4v 2 sin f3 2 
^m Ul m U2 



(Mi -M 2 )cos20 r 
■(Mi + M 2 )coBh26i, 



2v 2 sin [3 2 



(M l + M 2 ) cos0sinh26»i 



^/m Ul m V2 
2v 2 sin f3 2 



(Mi - M 2 ) sin0sin26» r 



(37) 



(3? 



We see that the imaginary part of the off-diagonal component depends both on the Majorana 
phase in U ([36|) and the phase 9i in R. Even if it could be measured, and the neutrino mass 
eigenvalues Mi j2 , m vi 2 were known, still only one combination of the angle factors 9 r ^ entering 
in leptogenesis (^) would be known, and there would still be an ambiguity associated with 
the Majorana phase 0. In fact, no CP violation is induced by the renormalization ([38|) in 
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this simple two-generation model, since it is not possible to define the Jarlskog invariant 
J I fl2~2|) and its analogues (|23|). Such invariants can be defined in a three-generation model, 



and CP- violating observables are demonstrably proportional to it, as we show in the next 
Section. 

As commented in subsections 2.3.3 and 2.3.4, non-degeneracy between the heavy singlet 

neutrinos iVj induces, via renormalization, dependences of the entries of m| and A e on phases 

in the product ZP 2 . In the two-generation case, this dependence is on the one phase in P 2 , 

since Z has no phases in this case. This makes changes in Arg(m|)i2 and Arg(v4 e )i2, but 

these are suppressed to the extent that log j^- <C log M ^ T ■ Moreover, these small changes 

are identical in the leading-logarithmic approximation. On the other hand, the corrections 

of the order of log j^- log M ^ T to the phases of the diagonal terms in A e may be sizeable. 

These are given in the two-generation model by 

r r + + ,ni x /m v ,m v Sm v , — m v ^)M\M 2 M> 
Im\[YjLY u ,YX] U ] = V 04 4o 2 log -/cosh 2^ sin 2g r sin 0, 

L J 2t> 4 sm p Mi 



^Jm Vl m V2 (m Vl + m I/2 )M 1 M 2 M 



+ — — — — " 1 " il — - — - log — - sinh 26 i cos 26L cos 6 

2t» 4 sin 4 /? & Mi y 

[YjLY u ,YX] 22 ] = -lm\[YjLY u ,YX] U ]- (39) 



Im 



We see that this effect vanishes if M 2 = M x . Defining M 2 = Mi (1+5), Eq.(p9|) grows with the 
dimensionless parameter 5 (linearly, if 5 is small) and is maximized when log[M2./Mi] = 1. 

3 CP- Violating Observables in the Charged-lepton 
Sector 

In this Section we discuss in more detail CP-violating and LFV observables in the charged- 
lepton sector. The slepton-mixing effects discussed in the previous Section generate LFV and 
CP-violating vertices involving charginos, which in turn induce effective non-renormalizable 
interactions, as we discuss in the following. 

3.1 Chargino and Neutralino Interactions 

The relevant neutralino and chargino interactions for leptons and sleptons are given by J?]. [33| 

£ = ^&x^ + ^P fl )x^ (40) 

where P R =(1 + 7 5 )/2, P L = (1 - 7 5 )/2 and 

N t L AX = -g{V2Un6 w (0 N ) Al (U e y xi+3 + ^ (0 N ) A3 (U e )* Xj }, (41) 

\J Im-w cos p 
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N&x = -9[-^={(0 N r A2 + tan9 w (0 N r A1 }(U e )* Xi 

+ J mlh R (o N y A3 (u e y w ], (42) 

\J Amy/ cos p 

C^ax = 9 J^ 13 J PcLYAUuYxj, (43) 
V 2mw cos p 

cf AX = -g(Oc R ) A1 (u„y Xi . (44) 

In these expressions the matrices On, Ocl, Ocr, U e and U u diagonalize the neutralino, left- 
and right- chargino, charged-slepton and sneutrino mass matrices, respectively. The indices 
run between A = 1, 4 for neutralinos, A = 1, 2 for charginos, X = 1, 6 for sleptons and 
X = 1, 2, 3 for sneutrinos. In our framework the complex phases appear only in U e and U u . 

3.2 LFV Muon Decays 



£ = ^{m^ARjI^a^eLF^ + m^AjjTZo^ e R F^ lv 



The effective Lagrangian for polarized /i + — > e + 7 and fi + — > e + e + e decays is [[33 

AG F 

_ _ 

+9i(m e L)(e R e L ) + g2{^Le R ){e L e R ) 

+fi'3(7AR7 / 'ei J )(e^7 M e i? ) + ^^T^X^Ta*^) 

+#5(7^7^) (el7/^L) + g&{Wl^^L){e^l^e R ) + h.c.}. (45) 

Here Al and A R are the dimensionless photon-penguin couplings which induce /i + — ► e^7 
and n + — > e R j, respectively, which also contribute to the /i + — > e + e + e~ process, and the 
gi : % = 1,...,6 are dimensionless four-fermion coupling constants which contribute only 
to /i + — > e + e + e~. Explicit expressions for A L>R and the g^ in terms of N^ AX , C^ AX are 



lengthy []33|, so we do not rewrite them here. 

In the notation (|45| ) the total fi + — > e + 7 branching ratio is given by 

5r(> + -> e+ 7 ) = 384tt 2 (\A l \ 2 + \A R \ 2 ) , (46) 

and that of fi + — > e + e + e~ by 

log-f- — \(C 5 + C 6 ) 
mj 4 J 

+16(C 7 + C 8 ) + 8(G ) + C 10 ). (47) 



The coefficients Cj appearing in fl4"7|) are functions of A L R and c^: 
16 16 
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C 3 = \g 5 \ 2 , C A = \g 6 \ 2 ,C 5 = \eA R \ 2 , C 6 = \eA L \ 2 , C 7 = Re(eA R g* 4 ), 

C 8 = Re(eA L g* 3 ), C 9 = Re(eA R g*), C 10 = Re(eA L g* 5 ) . (48) 

In order for CP violation to appear in any process, interference between different terms in the 
amplitude for the process must occur. Therefore, all possible observables in /i — > decays, 
such as differences between the fi + — > e + 7 and /i~ — > e _ 7 rates, vanish in the leading order of 
perturbation theory. Moreover, the process yT —>■ e~7 is not measurable with high accuracy 
because of the large backgrounds. However, when muons are polarized, a T-odd asymmetry 
for final-state particles in fi + — > e + e + e~ can be defined. Since CPT is conserved, the T-odd 
asymmetry measures the amount of CP violation in our model. 

The muon polarization vector P can be defined in the coordinate system in which the z 
axis is taken to be the direction of the electron momentum, the x axis the direction of the most 
energetic positron momentum, and the (z x x) plane defines the decay plane perpendicular 
to the y axis. It is necessary to introduce an energy cutoff for the more energetic positron: 
Ei < (m M /2)(l — 5). We use 5 = 0.02 to optimize the T-odd asymmetry, following [j33| . 
Assuming 100%-polarized muons the T-odd asymmetry is then defined by 

N(P y >0)-N(P y <0) 3 
At ~ N(P y > 0) + N(P y < 0) " 2Br(5 = 0.02) {2 -° Cu " h6Cl2} ' (49) 

where N(Pi > (<)0) denotes the number of events with a positive (negative) Pi component 
for the muon polarization, 

C n = lm(eA R gl + eA L g* 3 ), C 12 = lm(eA R g* 6 + eA L g* 5 ) , (50) 

and the branching ratio for 5 = 0.02 is 

Br{5 = 0.02) = 1.8(Ci + C 2 ) + 0.96(C 3 + C 4 ) + 88(C 5 + C 6 ) 

+14(C 7 + C 8 ) + 8(C 9 + C 10 ). (51) 



It is known [33] that the asymmetry (EH) may be large in SU(5) SUSY GUTs [34]. We study 



below whether this is also the case in the minimal supersymmetric seesaw model. 
3.3 Electric Dipole Moments 

The electric dipole moment of a generic lepton i is defined as the coefficient dg of the inter- 
action 

C = --d t to^tF^. (52) 
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The current experimental bounds are d e < 4.3 x 10 27 e cm for the electron ||35|| , d^ = 
(3.7 ± 3.4) x 1(T 19 e cm for the muon §6|, and \d T \ < 3.1 x 1(T 16 e cm for the r [|37|. 
An experiment has been proposed at BNL that could improve the sensitivity to d^ down 
10 -24 e cm fl3"8"|, and PRISM and neutrino factory experiments aim at sensitivities 



to d^ ~ 



dfj, ~ 5 x 10 e cm. These bounds will impose serious constraints on CP violation in the 
MSSM, as will prospective improvements in the sensitivity to d e and d T . 

In the MSSM, the di receive contributions from chargino and neutralino loops: 



d, 



+ df 



(53) 



where |L9, gO 



df 



(4vr) 2 



(4vr) 2 



23 m + ,m\ 

A=l X=l fx v x 

a-l A-l ly ly 



and the loop functions are given by 

A{r) = 
B(r) = 



1 



2(1 -r) 2 
1 

2(r- l) 2 



3 — r 



1 + r + 



2 logr 
1 — r 
2r log r 
1 — r 



(54) 
(55) 



where the relevant chargino and neutralino couplings were given above. We do not consider 
in our analysis the possibility of CP violation in the chargino and neutralino mass matrices. 



3.4 Role of the Jarlskog Invariants 

We first present approximate formulae for the effective couplings in (^5]), in order to show 
the qualitative behaviours of the LFV processes and demonstrate the role of the Jarlskog 
invariants. Since (5m|) and (5A e ) are the only sources of off-diagonal components, the 
only non- negligible terms are An, g± and g§: other terms are suppressed by the electron 
or muon masses. For illustrative purposes in this subsection only, we assign to the soft 
supersymmetry-breaking parameters common value irts at the electroweak scale: 

M 2 = M 1 = fi= (A e ) 22 /Y e2 = m s , 

{m\)ii = (m|)« = m| . (56) 
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Assuming mj ^> mz, we then find 

^ = 4.8xlO- 5 (^) 2 
x {Afi - O^A^A^ + 0.66A& - + tan/3(2.4Af 1 - 1.12A^A| l )} , (57) 

£ 4 = 6.4 x l0- 5 (l°g^) 2 
x {(1- 0.24 sin 2/5) A^ + (-0.81 + 0.12(sin2/3 + cos 2/3))Af3Af 1 } , (58) 

^ 6 = -1.9xlO- 5 (S) 2 
x |(1 — 0.70 sin 2/3) A^ + (—0.43 + 0.35(sin2/3 + cos2^))A| 3 A| 1 | , (59) 



where 



A' 



A e _ / (fiA e )ij/Y ei \ 



13 \ m s J 

We remind that A|i = (Af 2 )*. The sin 2(3 and cos 2/3 dependences of g 4 and g 6 above are 
due to Z-penguin diagrams. In the branching ratio for fi — > 3e, the contribution from Ar 
tends to dominate due to the phase-space integral. Then, assuming that Ar dominates in 
Br(ji + — > e + e + e~), the T-odd asymmetry At is given by 



An 



Im 


A^ A^ A^ " 

^12^23^31 


0.039 + 0.196 tan /3 + 0.017/ tan /3 


1 Al2l 


(1 + 2.4 tan/?) ^1(0.64 + 1.12 tan /3) 2 

^21 



(60) 



where we have expanded sin 2/? and cos 2/? in terms of tan/3. Also, in writing (|60|) , we have 
taken A^f = A^ 3 e = 0, for simplicity. 

We see explicitly how A T (|60f) depends on the Jarlskog invariant fl22"D , and it is apparent 
how analogous invariants Jai 2 > etc, with one or more Afj — > A^ e could also contribute. We 
see that A T could in principle reach ~ 10%. However, if ImfA^A^A^] <C |Af 2 | 2 , as 
one might expect, or if tan/3 3> 1, A T is suppressed. However, we stress that Eq. (j60|) is 
approximately correct only for (^) and cannot be used to predict At in more general cases 
which will be considered in the next Section. 

Next, we present approximate formulae for the effective coupling in ((5*2"|) in the specific 



case fl5*6"P. Since relative phases between (m|), (m|) and (A e ) contribute to the electric dipole 
moments, non- vanishing contributions to d\ come only from slepton diagrams, not sneutrino 
diagrams. From the explicit formula (1541), it is also clear that the sneutrino diagrams do not 
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give a non- vanishing value, since they depend on (U u )xi(U u )xi, and the CP-violating phases 
exist only in the mixing matrices of the sleptons, in our approximation. 

Since it depends on the neutrino model which contribution is dominant in the electric 
dipole moments, as shown in the previous Section, we first show the general formula for df° 
in the limit Q5fi|): 



e 



9l 



(An) 2 ms 



Im 



cjvAjjj A ili2 • • • A, 



N=l h 



■,IN 



(61) 



Here, Ay is the flavor-dependent part of the slepton mass matrix, which normalized by m|, 
and includes parts generated by the renormalization as well as tree-level parts. At least one 
of the Ajj in each product term involves the left-right mixing of a slepton. The coefficients 
d are: 



Cl 



1 

12'' 



(-2 



20 : 



C3 



30 : 



1 

42 : 



1 

56' 



(62) 



When the heavy singlet-neutrino masses are not almost degenerate, and corrections to the 
soft supersymmetry-breaking terms and the relative phases among (ml) and (A e ) are gen- 
erated at 0(log 2 M m n t )-> the approximate formula becomes 



9 2 v 



mi 



(47T) 



mt 



^tan/3 2 + ^tan/3 



105/ 



3 m 2 s jl 



30 (4tt)2 ml £ 



mt 



(63) 



where the first term is proportional to (|30"D and the second to (p9 ). 

When the heavy singlet-neutrino masses are almost degenerate and the correction to 
ImL4 e ]jj, proportional to log jj^ 1 log M ^ T , is dominant in the electric dipole moments, the 
approximate formula is 



9 2 y 



m t 



12 (4tt) 



Im 



(64) 



4 Numerical Analysis 
4.1 Calculational Procedure 



We first fix the gauge couplings, charged-lepton and quark Yukawa couplings and tan/3 
(which is a free parameter) at the scale Mz, and then run them with the two-loop MSSM 
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renormalization-group equations up to the scale Mjv- At M^, we introduce the heavy sin- 
glet neutrinos, fixing their masses, the light-neutrino masses and mixings according to the 
oscillation data. We then choose the matrix R and calculate Y v according to (0). Sub- 
sequently, we run all the Yukawa-coupling matrices from M N to M GUT using the one-loop 
renormalization-group equations M. At Mqut we assume universal boundary conditions 
( Effi) for the soft supersymmetry-breaking terms. We then run all the soft supersymmetry- 
breaking masses and Yukawa matrices back to M^, where the heavy singlet neutrinos and 
sneutrinos decouple. The soft supersymmetry-breaking mass matrices at low energies are 
obtained finally by running all the MSSM parameters back down to Mz- We use the elec- 
troweak symmetry-breaking conditions to fix the magnitude of the Higgs mixing parameter 
/i, taking its sign positive as motivated by — This sign is also consistent with the bounds 
from b — > sj. Then we calculate the squark, slepton, chargino and neutralino mass matrices 
and finally LFV rates, the T-odd asymmetry At in polarized fi + —>■ e + e + e~ decays and the 
EDMs of the electron and muon, for chosen values of the input parameters. 



4.2 Illustrative Results 



In our numerical examples, we take [39] Am 2 tm = 3.5 x eV 2 with sin 823 = 0.7 for 
atmospheric neutrinos, and Am^, = 5.0 x 10~ 5 eV 2 with tan 2 9 12 = 0.45 for solar neutrinos, 
corresponding to the LMA solution. Also, we fix tan 2 9 13 = 0.055, which is the largest value 
allowed by the global three- neutrino data fit. We parametrize the orthogonal matrix R in 
the form 



R 



' C2C3 — C1S3 — S1S2C3 S1S3 — C1S2C3 
c 2 S 3 C1C3 - hkh -S1C3 - C1S2S3 I , (65) 

\ ^2 S\C2 C\C2 



where #1, 9% , #3 are arbitrary complex angles. Using a generic R as input is crucial for 
neutrino phenomenology. We stress the following: 

(i) Because R is a complex orthogonal matrix, the values of its entries are not restricted 
to any small range, but rather are exponential functions of complex numbers. This 
implies via (|14"D that, for a suitable choice of #1,2,3, all the elements of Y v can be large 
even if one starts with small (for example hierarchical) neutrino masses; 

(ii) Large imaginary components are in general present in every entry of Y u . Since YjY u 
depends on the phases of both R and U, as seen in (|16D, sizeable CP- violating effects 
may be induced. However, if R = 1, only S in U contributes to YjY u . 
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We first calculate d e and using ( pi]) and (0), and study their possible ranges in our 
model by scanning over the allowed values of the free parameters. As expected, d e and 
d^ are very small. For degenerate right-handed neutrinos, typically d^ does not exceed 
10~ 31 e cm for the values of the parameters consistent with the bounds on fi — > cy, and 
is therefore many orders of magnitude below the sensitivity of any planned experiment. 
The approximate relation between muon and electron electric dipole moments in our model, 
dfj,/d e w — m M /m e , holds very well in this case. However, in the case of non-degenerate right- 
handed neutrinos, the \og(M N ./M N .) effects introduce a dependence of the electric dipole 
moments on the leptogenesis phases, as discussed previously. These new contributions may 
change the size of electric dipole moments by a few orders of magnitude. In particular, 



similarly to other supersymmetric models [40|, the sign of d^/d e may be altered and the 



naive relation d^/d e ph —m^/m e may be violated by a large factor. Detailed analyses of 
the electric dipole moments in the case of non-degenerate heavy neutrinos will be presented 



elsewhere 41 



On the other hand, the situation with the T-odd asymmetry in polarized /i + — > e + e + e~ 
decays can be different. As explained in the previous Section, the decay fi + — > e + e + e~ 
receives contributions from box diagrams and photonic penguin diagrams, with the latter 
usually dominating. However, if there are cancellations in the dipole-moment-type fi — e — 7 
vertex, the box and penguin contributions to fi + — > e + e + e~ may become comparable. In 
that case, if there are large CP-violating phases present in the slepton mass matrices, the 
T-odd asymmetry At can be large. This implies an anti-correlation between Br(fi — ► cy) 
and At ■ the latter can be large only if the former is suppressed. 

To illustrate how such a cancellation might come about, consider the two-generation 
example described at the end of Section 2, which also applies in the full three-generation 
case if m U3 <C m Ul2 : we assume for simplicity that also m U2 <C m Ul . We see from (|37D that 
(YjY^j is suppressed if 



(Mi + M 2 ) cos 26 r w (Mi - M 2 )cosh2^, (66) 
and that the smaller quantity Im (YjY v ^ is also suppressed if 

(M l + M 2 ) cos 0sinh2#i w (M x - M 2 ) sin sin 26 r . (67) 



Both the conditions ( p6lp7|) may in principle be satisfied simultaneously for suitably tuned 
values of 9{, 9 r and (f>. 

In a general and more physical case, the cancellation in Br(fi — > erf) is much more 
complicated and depends on all free parameters, namely the phases and mixing angles in 
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fl — > e e e 




Figure 1: Branching ratios for the decays (a) /i + — ► e + 7 and (b) fi + ^ e + e + e~ as functions 
of the common soft mass mo and the Majorana phase <p2, for the fixed choice of neutrino 
parameters described in the text. 

the matrices U and R, the values of the light- and heavy- neutrino masses, the choice of the 
soft supersymmetry-breaking initial conditions, details of the renormalization-group running 
procedure, etc.. However, the fact that such cancellations occur is robust and qualitatively 
well understood. To give a representative numerical example, we choose hierarchical neutrino 
masses with m Ul = 0.028 eV, M Nl = 1.2 x 10 15 GeV, M N , 2 = 1.5 x 10 15 GeV, M Nz = 3 x 10 14 
GeV, 5 = 7r/2, 0i = —1, 6*i = 0.3z, 9 2 = 0.5z, 9 3 = O.li, and use the same neutrino 
oscillation parameters as above. We run the renormalization-group equations with right- 
handed neutrinos down to the scale = 3 x 10 14 GeV. 

Fixing these parameters, we first choose mi/2 = 500 GeV, Aq = GeV and tan/5 = 20 
and scan over the remaining free parameters mo and <p2- In Fig. [I] we plot the branching 
ratios of the decays [i + — > e + 7 and [i + — > e + e + e~ as functions of the common sfermion soft 
mass parameter m and the Majorana phase 02- For a large region in the plotted parameter 
space, Br(fi + — > e + 7) is below the present experimental bound, whilst Br(fi + — > e + e + e~) 
does not at present impose any constraints. The branching ratios in Fig. [I] are correlated, 
implying that the photonic penguin diagrams are dominating also in fi + — ► e + e + e~ . For the 
plotted values, both rare decays should be observed in the planned experiments. 
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T-odd asymmetry A T 




Figure 2: The T-odd asymmetry A T in polarized jj, + — > e + e + e decay for the same set of 
parameters as in Fig. 

The T-odd asymmetry At is plotted in Fig. [| for the same set of parameters as in Fig. |l[ 
Comparison with Fig. [I] shows that At is strongly anti-correlated with the branching ratios. 
For the region with the deepest cancellation, the T-odd asymmetry is negative. Whilst its 
absolute value may exceed 10%, the branching ratio of fi + — > e + e + e~ is relatively small 
in that region, making its precise determination difficult. However, for large positive At, 
Br(fi + — > e + e + e~) exceeds the 10 -14 level, implying that several hundred events could be 
observed in the planned experiments. Since the experimental sensitivity to At should scale 
as 1/ a/ N events , these future experiments would be able to measure a non-zero value of At 
for large ranges of mo and 02 in Fig. |2|. 

Whilst the long-baseline oscillation experiments at neutrino factories will measure the 
phase 5 in (p]) , we stress here again that the T-odd asymmetry At depends on all the 
phases in the matrices U and R. Thus different combinations of phases in the Yukawa 
matrix Y v are probed in the neutrino oscillation and stopped muon experiments. This point 
is made explicitly by the dependence of At on the Majorana phase 02 in Fig. 0. 
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Branching ratios 



T-odd asymmetry A- 




Figure 3: (a) Branching ratios for the decays fi + — ► e + 7 and fi + — > e + e + e~ and (b) the 
T-odd asymmetry At in fi + — > e + e + e~ decay, as functions of the common soft mass mo, for 
the fixed choice of neutrino parameters described in the text. 



We have chosen the neutrino parameters in such a way that the absolute values of all 
the neutrino Yukawa couplings are close to unity. This induces large rates of LFV and CP 
violation. However, due to the cancellations in the photonic penguin diagrams, the induced 
Br(fi + — > e + 7) can be consistent with the current experimental bounds even for relatively 
small sparticle masses. In Fig. ^|we plot (a) Br(fi + — > c +r )) and Br({i + — > e + e + e~) and (b) 
At as functions of mo for fixed mi/2 = 200 GeV, Aq = GeV, tan/3 = 10 and 02 = 2.1. For 
a small region around mo = 300 GeV, Br(fi + —>■ e +r y) is below the present limit. At the same 
time, the T-odd asymmetry in [i + — > e + e + e~ may be as large as 10% for the allowed values 
of mo- The dependence of the branching ratios and the T-odd asymmetry on the Majorana 
phase 2 is demonstrated in Fig. |j. Again, large A T is expected if the decay /i + — ► e + 7 is 
suppressed due to the cancellation. 

Finally, we note that the branching ratio of the decay r + — ► /i + 7 does not have cancel- 
lations in the parameter region considered, as seen in Fig. [5[ Therefore, in this scenario 
Br(r + — > /U + 7) might be just below the present experimental bound and discoverable at 
the LHC or the B factories. As the decays r — > Hi are suppressed relative to r + — ► /U + 7, 
detailed quantitative studies of them are beyond the interest of the present work. 
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T-odd asymmetry A- 




Figure 4: (a) Branching ratios of the decays fi + — ► e + 7 and fi + — > e + e + e~ one? f&j £/ie T-odd 
asymmetry At in fi + —>■ e + e + e~ decay, as functions of the Majorana phase <p2 for itlq = 300 
GeV. All other parameters are fixed as in Fig. 



5 Discussion and Conclusions 

We have seen in this paper that the T-odd asymmetry At in polarized fi + — ► e + e + e~ 
decay may offer the best prospects for studying in the laboratory CP-violating effects in 
the minimal supersymmetric seesaw model. On the other hand, the electric dipole moments 
of the electron and muon are suppressed in the minimal supersymmetric seesaw scenario 
discussed here. This is because the CP- violating phases are induced by renormalization-group 
running only in the off-diagonal entries in the slepton mass matrices, as already discussed 
The naive relation d^/de ~ —m^/m e holds very well in the case of degenerate 



in 



right-handed neutrinos. In the case of non-degenerate right-handed neutrinos, logarithmic 
effects arising from \og(M N ./M N .) introduce a dependence on the leptogenesis phases. These 
new contributions may become dominant and the naive relation d^/d e ~ —m^/m e is badly 
violated. In the most optimistic case, the electric dipole moments may approach the level 
observable at the proposed experiments ET . 



The possibility of measuring a non-zero At could have far-reaching consequences, since 
it provides complementary information on the CP-violating phases in the neutrino Yukawa 
matrix Y v . As has been discussed, At depends in the leading- logarithmic approximation on 



26 



Branching ratio 



u K10 

ffl 0.59 
0.58 
0.57 
0.56 
0.55 
0.54 
0.53 
0.52 
0.51 

280 290 300 310 320 330 

m [GeV] 

Figure 5: Branching ratio of the decay r + — ► fi +/ y as a function of mo. All other parameters 
are fixed as in Fig. [|. 

a single combination of light-neutrino phases and the three phases in YJY V that contribute 
to leptogenesis, whereas the CP-violating phases in the light-neutrino effective mass matrix 
depend on the phases in Y„Yj, that do not contribute to leptogenesis. 

This is one reason why At may be observable even if CP violation is undetectable in 
neutrino oscillations. We recall also that the latter is in practice observable only if the 
neutrino masses and mixing angles are favourable. For example, if ~ and/or Am^, 
is small, as in the case of vacuum oscillations, CP violation is unobservable using neutrino 
factories. However, as seen in (|I^), (some of) the Yukawa couplings in Y v may still be large 
and imaginary, implying that At might be large. 

On the other hand, a large value of At requires cancellations in the slepton-induced 
fi — e — 7* vertex, which happens only in a restricted region of the parameter space. The 
asymmetry At is anti-correlated with the branching ratio of fi — > ej, and it can reach ~ 10% 
if fi — > cy is suppressed. The asymmetry At may be measurable in planned high-intensity 
stopped-muon experiments, which aim at a sensitivity to Br(/i — *> eee) ~ 10~ 16 . 

In the case of r — > /ry the cancellation does not happen for the same parameters as in 
\l — > erf. Therefore, in the scenario considered in this paper Br{r — > ^7) is large and can be 
observed at the LHC experiments. 

It is interesting to review what we would learn if non- vanishing At were observed. The T- 
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odd asymmetry At is approximately proportional to the CP invariant (|22|). We recall that, 
if Y v has a hierarchical structure, the off-diagonal components of the left-handed slepton 
mass matrix are given by 

H) 12 « (y°) 2 x 31 x; 2 + (y u d Jx 21 x* 22i 

The Jarlskog invariant may have a sizeable value if {Y^) 2 X^\X^ 2 and (Y®) 2 X2iX 22 are 
comparable. Thus, if non-vanishing At is observed, the generation structure in the neutrino 
Yukawa coupling may be constrained, as well as the CP-violating phase. 

In conclusion: searching for CP violation in lepton-flavour-violating processes is a pos- 
sibility that should not be neglected, since it provides information complementary to that 
provided by neutrino oscillation experiments. In particular, At may be measurable even if 
CP violation is unobservable in neutrino oscillations. 
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